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ABSTRACT 


Addition of the direct and indirect planetary perturbations to the 
author’s 1963 modification of Hansen’s theory is proposed. The method 
described can determine the Saturnian effects upon the motion of the 
outer Jovian satellites. The expansion of the disturbing function and of 
its derivatives is reduced to a form convenient for programming. 



BASIC NOTATIONS 


the jovicentric position vector of the sun 

the jovicentric position vector of the fictitious sun moving in Hansen T s mean 
ellipse 

Hansen’s perturbations of the solar radius vector 

the elevation of the sun relative to the plane of the Hansen mean ellipse 
Gibbsian vectors of the solar mean ellipse 

the semimajor axis, eccentricity, and mean daily motion associated with 
the solar mean ellipse 

the undisturbed mean anomaly of the sun 

Hansen’s perturbations of the mean anomaly of the sun 
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refer to the heliocentric motion of Saturn; and the notations 
r, r, v, u, P, Q, R, a Q , e Q , n Q , 

g " n o t + So’ n o Sz ’ ^ = g + n 0 7 

refer to the jovicentric motion of the satellite 
the mean argument of the perigee of the satellite 

the mean longitude of the ascending node of the satellite. The plane of the solar 
mean ellipse is taken as the basic reference plane. 

the mean inclination of the orbital plane of the satellite 
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A'.B'.C 
A" ,B",C 
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on P, 


. ] - the matrix which is a polynomial in Euler parameters k ± , k 2 , \ 3 , k and which 
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' - the projections of r" 

: - the projections of D 
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PLANETARY EFFECTS IN THE MOTION 
OF NATURAL SATELLITES 

by 

Peter Musen 

Goddard Space Flight Center 


INTRODUCTION 

In this article we develop the theory of general planetary perturbations in the motions of the 
"natural" planetary satellites— one which can be applied, for example, to the outer Jovian satel- 
lites to determine the effects of Saturn. This topic represents a particular case of the four-body 
problem. 

For the determination of the purely solar effects in the motion of the outer Jovian satellites, 
the author has suggested a modified form of Hansen’s theory (Musen, 1963). In the motion of these 
satellites, the solar effects are strongly dominant over the effects of Saturn or any other additional 
small force. Thus the determination of the direct solar effects, under the assumption that the 
solar motion is Keplerian, constitutes the "main problem" of the theory of satellite motions. 

A modified form of Hansen’s theory as suggested by the author was programmed by Charnow 
(1966), and as a test the actual expansion of the solar perturbations and the ephemerides for 1967 
and 1968 were computed for Jupiter X. The satellite was found by E. Roemer in close proximity to 
the predicted position. The representation, by this theory, of 30 years (1968) of Jupiter X observa- 
tions suggests that it may be necessary to include Saturnian effects in the theory of motion. 


We must expect the terms affected by small divisors in the integration process to be among 
the most significant terms in the expansion of the planetary perturbations. Such terms, in the 
past, were discovered through attempts to find critical linear combinations of the mean anomalies 
of all three bodies (i.e., combinations with very small mean motions). Obviously, such an approach 
requires luck, and reasonable doubts can be raised that all the significant planetary effects can 
be found in this way. We prefer to use a form of the theory and a computational scheme which 
permit the computer to select the significant terms automatically and to decide their importance 
on the basis of a purely numerical criterion. 

Two types of planetary effects are to be included in the theory: (1) the direct planetary 
effects caused by the direct perturbative action of Saturn on the satellite, and (2) the indirect ef- 
fects caused by the deviation of the solar motion from its Keplerian approximation. The indirect 
effects must be included in the solar part of the disturbing function. Thus they represent an im- 
mediate contribution to the main problem. In order to include the direct planetary effects in the 
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theory presented here, necessary changes are made in the formulation of the main problem. The 
basic ft small parameter" in front of the solar part of the disturbing function is larger than that for 
the planetary part. This means that the influence of the mutual perturbations of Jupiter and Saturn 
is greater in the indirect than in the direct planetary perturbations of the satellite. Under the in- 
fluence of the planetary effects, the mean argument of the perigee and the mean longitude of the 
ascending node of the satellite will cease to be linear functions of time. They will also contain 
terms that are quadratic in time; and, theoretically, also terms of higher orders. Thus, under the 
influence of the planetary effects, the motion of these two angles will be slightly accelerated. 

The form of the theory we are pursuing here is numerical. The values of the mean elements 
of all bodies can be substituted from the outset, and the final output will be trigonometric expan- 
sions with purely numerical coefficients. We employ ideas expressed by the author in his earlier 
work on the lunar perturbations of artificial satellites (Musen, 1961) and expand the disturbing 
functions into a series of polynomials in A, B, C, where A is the projection of D, the jovicentric 
position vector of Saturn, on the mean line of apsides of the Jovian satellite; B is the projection of 
D on the direction normal to the mean line of apsides in the orbital plane; and C is the projection 
of D on the line normal to the satellite’s orbital plane. These polynomials are very simple and 
their trigonometric expansions can be obtained fairly easily with a computer. 

The disturbing function for the direct planetary effects contains the odd negative powers of f D I . 
Following Newcomb’s idea in his work (Newcomb, 1907) on the planetary inequalities in the motion of 
the moon, we abandon the application of Laplace coefficients to obtain the trigonometric expansions 
of the powers of | D| , and suggest instead the use of harmonic analysis. Of course, the Keplerian 
elliptic approximations to the jovicentric motion of the sun and to the heliocentric motion of Saturn 
will produce the most significant part in the expression for the direct planetary effects. However, 
we shall provide a device to carry the process to higher approximations, either in order to include 
them, if necessary, or to have a more precise idea of the magnitudes of the terms we omit. In this 
work, in accordance with Hansen, we use the satellite’s osculating orbit plane as the basic reference 
plane for the expansion of the perturbations in the satellite’s radius vector and mean anomaly. In 
the classical Hansen theory, the motion of the satellite’s osculating orbital plane is referred to the 
osculating solar orbital plane. However, the elevations of the sun and of Saturn relative to their 
mean orbital planes are very small, so that for our purposes we can refer the motions of the 
osculating orbital planes of the sun to the mean orbital plane. Also in this work we can use the 
existing Hill’s expansions (Hill, 1890, 1906) of the Hansen coordinates of Jupiter and Saturn. 

THE DIRECT PLANETARY EFFECTS IN THE ORBITAL PLANE OF THE SATELLITE 

We take the mean orbital plane of the sun relative to Jupiter as the basic reference plane. The 
x-axis is directed toward the perijove of Hansen’s mean ellipse of the sun, and the z-axis is nor- 
mal to this plane. We use the following notations: 

r - the jovicentric position vector of the sun, 

r' - the jovicentric position vector of the fictitious auxiliary sun moving in 
Hansen’s mean ellipse, 
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V 1 


u ' 


P\ Q\ R' 


t + 


n o # Sz ' 


l' - g' + n Q f Sz' 


f ' 


Hansen’s perturbations of the solar radius vector, 

the elevation of the sun relative to the plane of the Hansen mean ellipse, 
Gibbsian vectors of the solar mean ellipse, 

the semimajor axis, eccentricity, and mean daily motion associated 
with the solar mean ellipse, 

the undisturbed mean anomaly of the fictitious sun, 

Hansen’s perturbations of the mean anomaly of the sun, 

Hansen’s disturbed mean anomaly of the sun, 

Hansen’s disturbed true anomaly of the sun. 


The notations 


r " , 


r , 


i '“ , u" , P , Q , R' , 


g" 


g" * n 0 ' Sz ’ f ' 


refer to the heliocentric motion of Saturn; and the notations 


r, r, i/, u, P, Q, R, a Q , e Q , n Q 

g ‘ n 0 t * g 0 , n 0 ,v z ■ l ‘ s fr V Sz - f 

refer to the jovicentric motion of the satellite. The basic relations between the two types of posi- 
tion vectors (for example, r" and r") are: 


r" 

(1 + v" ) (r " + R" u" ) , 

(i) 

r " 

= P" r" cos f" + Q" r" sin f" , 

(2) 

r ' 

= (1 + v' ) (r 1 + R' u' ) , 

(3) 

r 

= P' r' cos f ' + Q' r' sin f ; , 

(4) 

r 

= ( 1 + v) r , 

(5) 

r 

= Pr cos f + Q r sin f . 

(6) 
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The expansions of P, Q, R into periodic series can be obtained from the relation 


[P, Q, R] = A 3 (-*/) • A ■ A 3 (+w) (7) 

given by the author in his previous work (Musen, 1963), where A 3 (a) is the matrix of rotation 
around the z-axis; the A-matrix is a polynomial in the Euler parameters \ , A 2 , A 3 , A 4 which carry 
all the periodic oscillations of the osculating orbit plane around its mean position; andw, -a>' are 
the mean argument of the perijove and the mean longitude of the ascending node respectively. 

We obtain from Equation 7: 


P x - + _ ^^2 ) cos ) _ 2A 3 A 4 sin (a, - f (\j 2 - A 2 2 ) cos (a. + a;' •) - 2Aj Aj sin (to- + a ' ) , (g) 

P y = + (\ 4 2 - A 3 2 ) sin (« -w' ) + 2A 3 A 4 cos (a, - a;' ) - - \ 2 2 ) sin (oj+ - 2Aj A 2 cos (a; + a;' ) , (9) 

P 2 = + 2(A 2 A 4 ■* Aj A 3 ) cosw + 2(Aj A 4 - A 2 \ 3 ) sina. , ( 10 ) 

Q X - - (\ 4 2 - \ 3 2 ) sin (., -a.' ) - 2 \ 3 V 4 cos (a - a' ) - \ sin (gj + u ' ) - 2 -Vj A 2 cos (a> + a' ) , (H) 

Q y ■ ■> (a 4 2 - V 3 2 ) COS (u - 2\ 3 V 4 sin(a -a ')- (Aj 2 -A 2 2 ) cos (a- a') - 2 A.J A 2 Sin (Ct a') , ( 12 ) 



+ 2(\ t 

^4 

■ss) 

COS 6J 

~ 2 (\ ^3 + ^2 \) sina; ’ 

(13) 

R 

X 

+ 2 (\ 

S 

■d- 

CM 

1 

cosw' 

- 2(\j A 4 4- A 2 A 3 ) sin a-' , 

(14) 

R 

y 

- 2 (S 

V 4 

*ss) 

COS (4 ' 

“ 2 (S V 3 - V 2 A 4 ) sine' . 

(15) 


R 2 = - \ 2 - K 2 + ^ + V • (16) 

We can set 



• 0 + XX 

sin -r£ + oK 1 , 

(17) 

S 

r S S • 

(18) 

^3 

= &a 3 . 

(19) 

\ 

Io + S X 

cos ~2 + o \ 4 , 

(20) 
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where i 0 is the mean inclination of the satellite's orbital plane and S\ 2 , 8\ 3 , S\ 4 are of the 
order of perturbations. Making use of Equations 17 through 20, we can write Equations 8 through 
16 with sufficient accuracy in the form 


+ S«/r x 

P o • 

(21) 

+ X 

0o . 

(22) 

+ S</i x 

R 0 , 

(23) 


where 


T o ! 0 

+ cos 2 ~2 cos (to - to' ) + sin 2 -j* cos(6j + o/) 


+ cos J 


2 sin (r/j Ctj ) 


sin 2 ~2 sin (&.' + &’') 


sin I sin co 


(24) 


r 


and 


sin ( , - ' ) 


~ u , 

m 2 s in (oj + ^ ) 


+ cos'* 


u . _ u 

~2 CO s ( ca ~ co ) - sin 2 -?£ cos (oj + a; ) 


sin I Q cos co 


J 


Si/r 


(25) 


sin I Q sin w 1 
sin I Q cos u 1 > 

cos i 0 J 

are the main long-period parts in p, 0, and R, and 

(cos ~2 SX _ 1 - sin -y S \ 4 ) ( P Q cos co - Q q sin co ) 

- (cos ~2 SX 2 - sin 8 \ 3 ) ( P Q sin co + Q q cos co ) 

( h r o \ 

+ ^ sin ~2 S \ 2 + cos ~2 S^ 3 jR 0 


(26) 


= P 0 8 *1 + 0o ^2 + R 0 S ^3 • 


(27) 
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Designating by D the jovicentric position vector of Saturn, we have 


D - r " + r ' ; 


( 28 ) 


or, taking Equations 1 and 3 into account, 


D - (l + v")(r"+R"u") + (l+v')(r'+R'u'). (29) 

In order to achieve symmetry of formulas and computational procedures, we decomposer", r',and 
D along P, Q, and R : 


r" = A" P + B" Q + C" R , (30) 

r' = A' P + B' Q + C' R , (31) 

D = AP + BQ + CR , (32) 

A = A" + A' , B - B" + B' , C = C" + C' . (33) 


The disturbing function H" associated with the direct Saturnian action, when expanded into a series 
in Legendre polynomials, can be written in the form 



where m" is the ratio of the mass of Saturn to the mass of Jupiter. Making use of Equations 32, 
and substituting 5 and 6 in the last equation, we obtain 
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where now m" designates the ratio of the mass of Saturn to the mass of the sun. 

A similar form of the expansion of the disturbing function was used by the author (Musen, 1961) 
in the theory of the lunar effects in the motions of artificial satellites. The application of harmonic 
analysis in g' and g" would be the easiest way to expand D -3 , D~ 5 , D~ 7 , • • • into trigonometric series 
In order to make this harmonic analysis possible, we must expand a Q 0" in powers of v' 9 v" , u', u" 
and Si/r. 

Setting 


D 


o 


+ r 


(36) 


5D - z>" r " + v ' r ' + R" u" + R' u ' , 


(37) 


we have 


D - D 0 + SD + H ^ D 0 , (38) 

where D 0 and also and SD are referred to the frame (P 0 , Q 0 , R 0 ). From now on there is no need 
for the subscripts on P Q , Q 0 , R 0 , and D 0 . We omit them and use simply P, 0 , R, and D. The notations 
A" , A ' , A ... will be associated with this new meaning of P, Q , R. We shall also change the 
meaning of 


a." sin * ’ „ i sin ^ 

0 a o 

and define them as the coordinates associated with the undisturbed Keplerian motions, depending 
upon g" and g' respectively. Thus 


r" _ 

cos f " - cos E" - e " , 

a o" 


r ' _ 

cos f ' - cos E' - e ' , 

a o' 


7 



sin f" 

a 0 

- e 0 " 2 sin E" , 

rr sin f ' - 

a o 

- e Q ' 2 sin E 

E" - e Q " sin E" 

= g" . 

E' - e Q ' sin E' 

= g' • 


Now, omitting the factor 1 + v and assigning the new meaning to A, B, c and D = r" + r' , we set 


a o 



1 

2 



2 


+ 


3 

4 



A 2 - B 2 



cos 2f + 


3 

2 





sin 2f 



A 2 + B 2 



3 

2 



cos f 


B /l5 
a o \ 8 D? 


A 2 + B 2 


3 a o" S \ / r \ 3 . 

o sin f 

^ D 5 / \ a o/ 


a o'' 7 a(a 2 - 3B 2 ) ( r \ 

8 D 7 ' a; 3 UJ C ° S 


B(3A 2 - B 2 ) 


sin 3f 


' (39) 


Designating the old value of a 0 f!" by a Q OJ, we have 


3 0 n o" - a 0 n " + K 3 0 n " 


(40) 


where 


D a o" d \ 

a o ' D dD J' 


and 


V = P + 0 — 7 + R —7-7 • 

D/a 0 » <9 A/ a Q v <9B/a Q <9C/a 0 " 


K = 


c , 

vr — - + n Sz 
t? r 0 


d „ , <9 

+ n n 8z + 

<?g 0 ^g" 


St/>> 


a„' 


d 4 „ 


(41) 


( 42 ) 
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Taking Equations 28 and 37 into account and neglecting the very small quantities u' 2 , u" 2 , u' 
u ' v ' , u" we have 


_ a d 6 

K = 'vx — =r + n ' Sz ' - — 7 + n " Sz" 




<9 g 


<9g" 


( a o / JD^ d \ 

a o" a o' V D/a 0 + a o" D dD/a o) 


„ r_^ ' I J)_ a _^ d 

a o" \ D/a o" + a o" D dWa o 


d o u' / r" ^ a 

+ a o" a o' \ n/a o + a o" ° 5D/a o" 


u 

+ — R 
a n 


V , + — — 

D/a" 


a ° r * _ 
C < 9 D/a„ 


80 • — • x V 

a o D/ «c 


We introduce the following auxiliary notations: 


P * P ' [P x ' cos (x ~ c c' ) + P y sin ( x - % 1 )] 


[ P x ' cos (oj + a' ) - P y " sin ( / + .*/ )] sin 2 + P ? sin I Q sin' 


JU ^ 

2 " = Q" • P - [Q x " cos (a;~ a;' ) + Q^" sin (oj- a)' )] COS 2 -y 


[Q x " cos (u> + oj' ) - Q y " sin (a. + u;' )J sin 2 + Q z “ sin I fl sin c 


a 3 “ - R" • P - [ R x " cos (w - oj' ) + R y " sin (6 j - o' )J cos 2 -j- 


+ [r; cos (tj + ) - R y " sin (m + to' )J sin 2 -j- + R z " sin I Q sin < 


(43) 


(44) 


(45) 


. (46) 
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b l " = P" • Q = [ P y " cos (a; - a/ ) - P x " sin (a; - a/ )] cos 2 - 7 ^ 


[ P y " cos (a ; + a/ ) + P x " sin (a; + co 1 ) j sin 2 + P z " sin I Q cos oo , (47) 


b 2 " = Q" • Q = [Q y " cos (cj - a>' ) - Q x " sin (a; - a>' ) J cos 2 


[Q y ' cos (a; + co' ) + Q x " sin (a; + a;')] sin 2 -y + Q z ' sin I Q cosw , (48) 


h> 3 " - R' • Q = [ R y ' cos (co - a/ ) - R x " sin (a, - co' )J cos 2 -j" 


[ R y " cos (a; + co' ) + R x " sin (co + a;' ) J sin 2 +■ R z " sin I Q cos co , (49) 


c " = P" * R - - ( P y " cos + P x " sin a/ ) sin I Q + P z " cos I Q 


c ^ = Q * R = - (Q " cosa;' + Q x " sin a' ) sin I Q + Q/ cos I Q 


3 " = R" * R = - ( R y " cos a:' +R x " sin a;') sin I Q + R z " cos I Q , 


= P' * P = P 


Q' • P = P 


R • P - P. 


" Q x 


COS 2 ~2 cos (co -Cl') + sin 2 ~2 cos (co + co') 


I 0 . 2 J o 

- cos 2 ~2 sin ( 0 . - 0 .') - sin 2 sin (a, + a ) 


sin I Q sina , 


cos 2 ~2 sin (a;- a;') - sin 2 sin(a;+a;') , 


b/ = P' ‘ 0 

b' = Q' * Q = Q y = + cos 2 cos (co - co ' ) - sin 2 cos (w + ) 

b 3 ' = R' • 0 = sin I Q cosa , 


l n I. 

T 2 


S' = p ' • R = R * 


sin I Q sin co' 


Q' • R - R y - - sin I Q cosa/ 


(50) 

(51) 

(52) 

(53) 

(54) 

(55) 

(56) 

(57) 

(58) 

(59) 

(60) 
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C 3 ' = R' H = cos I 0 . 


(61) 


It is of interest to note that these auxiliary quantities are affected by the long-period perturbations 
only. We have, also, 


A" 


B" 


a" cos f" + a " — - sin f " 

^ cx 


b," — - cos f" + b," — ~ sinf" 
a o a o 


(62) 

(63) 


C" 


c," cos f " + c " sin f“ 

a o a o 


(64) 


V 

a " 


o / r — r — 

i a/ — r cos f # + a ' — sinf' 


(65) 


B' a o l , V' , r' 

— — “ lb. — ; cos f + b_ — sinf 

a o" a o" \ a o a o 


(66) 


c: 


0 , r T, , 

C , — 7 cos f + c« 


r' . 7 ,\ 

— sinf ) 
a A / 


(67) 


A 

a o" 


A" + A' 

a o" 


B B" + B' 


C" + C # 


(33) 


After we have obtained the expansions in Equations 62 through 67 we deduce the expansion 


a o # r' D AA' + BB' + CC' 



( 68 ) 



AA" + BB" + CC" 


(69) 



P " 


Q " — - sin f " , 

Z a o 


(70) 


R" • 


r 

a ' 


(71) 
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All the expansions in Equations 53 through 71 are purely periodic, and depend upon the basic 
arguments g', g and a-'. 

The scalar expression of the operator K to be used in the actual computations takes the form 

_ a , o , <9 „ e „ d 

K - v r — rr + n bz ; + n n bz 

dr 0 <9g' 0 c9g" 


, . A' d B' (9 

+ v' I + 


a o JA/ a q a Q " <9B/a 0 " 


a o r ' D a o r) \ 

a o" a o' a o" D - 9D/a o" / 


A" <9 B" <9 r" D a o .9 


a 0 " <9A/a 0 " a Q " <9B/a 0 " a 0 " a 0 " D <9D/a 0 " 


3 0 U ' / <9 r? 

— — P + 0 

a 0 " a; V 2 dA/a Q Vz <9B/a Q 


n/ . ^ * 

B 7 n 


a 0 " D ^D/a 0 " 


u " I „ a 

a, — + b. 


a o" 


0 r' a 0 (9 

+ — R • — IT — 


a „• “ m/a" 


C <9 C <9 

by. + b'JJ 

a Q " 1 <5 B/a a 0 " 2 SA/a,," 


, . A <9 B <9 

— 


3 Va o" f7B a o" a o" ■ 9A a o" 


(72) 


The expansions of b +1 , 6u 2 , dv 3 to be substituted into Equation 72 can be obtained from Equation 27: 



( ^ 0 


V, \ / 


o 
t— 1 

fci/q - 

+ V cos T 


- sin ~2 ) cos ^ “ \ 

v cos ~2 

- sin ~2 bk 3 ) 


/ !o 



( „ 

l ° c ■ 

H 2 = 

- ^ cos ~2” 


- s in ~2 ~ 

^cos ~2 ok 2 

- sin — oK 


*3 = 


T U 2 + cos T iK 3 


(27') 

(27") 

(27’”) 


Not all the perturbations in the expansion of k (Equation 72) are comparable in size. For each 
satellite we shall make a special decision about their relative importance, and in each particular 
case the expansion of the operator K and the result of its application to a 0 fi" will always be 
reduced to only a few terms. 
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We obtained the differential equations governing the variations of Hansen’s elements h Q /h, 


T = 2 n 


e cos (x -77 0 ~ n o yt ) “ e o 


1 — 0 2 

1 e o 


'P = 2 TT 


h e sin (x~ rr o ~ n o V 4 ) 




by combining ideas from Hansen's planetary (1857-59) and lunar (1862) theories (Musen, 1963). We 
have, in the general case. 


dT 


* a o n 


19 a o n 


dF = +n 0 y * + M i + N i r ^97 


d<f 

dt 


n o y \ *~o r~, 


2 1 + M 2 


<7 a o n 


<9-5 


+ N 2 r 


S a 0 n 
<9 r 


_d_ ^0 
dt h 


N 3 r 


<7 a p b 
777 


where M t , N, (i = 1, 2, 3) are defined by the equations: 

'M'-i- 6) • (-E? - *) r. f. (' - d 


2 n „ 


2 n 


o r sin 


1— ‘‘o 2 ■'‘o/i-f'o 2 

rb K 2 ^ 


v a o /h 2 \ / a o \ 

7-b -Ob- 1 ) 


1 -e Q 2 


2n . 




h 2 ' - 2 

COS f • 3ejd 2 - 7 sin f ■ ( ~ ; - 1 


h o 2 


sin f 
a 2 1 “ e 2 


( 1 _e o 2 ) 3 


£ cos f - 2 e 0 + /l - e 0 2 ^ sin f 


til 

h o 2 


r sin f d 

a ° yi-e 0 2 d * a ° 


(73) 

(74) 

(75) 


(76) 


(77) 


(78) 


(79) 


1 - 0 2 2 * 

1 e 0 a Q 

(80) 

n o 7 e o sin7 


i ~ e o 2 a ° ,/! -e 0 2 

(81) 
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Now, however, in order to separate the main part of the direct planetary effects from a less sig- 
nificant part, it will be convenient to change the meaning of some symbols again. So far the nota- 
tions 7 and T have meant the coordinates which are functions of the disturbed mean anomaly 
Now we shall designate by the same symbols the coordinates which will be functions of the undis- 
turbed mean anomaly g. Thus 


— COS f 

a o 


cos E ~ e Q , 


~T~ sin f 
a 0 


yi - e Q 2 sin E ; 


and E-e 0 sinE - g. This change in notations requires the expansion of the right sides of Equa- 
tions 76 through 81 in powers of n Q bz. 

We now set 


2n 


l-e 0 2 e o 


0 1 ( o r ‘ 

— 1 -e 2 - 


(82) 


2n„ 


iM, 


1 - e Q 2 




(83) 


2n 0 r _ 
Nl = (1 -e Q 2) 3/2 ^ Sln f ’ 


(84) 


2n, 


iN, 


(1 


~ e Q 2 ) 3 2 




-L sin f 


(85) 


(l-e/) 3/2 .f( 2 ^ 3e 0 ) dg , 


( 86 ) 


2n n 




_ * 2 


1 - e 


r _ 7 sin f / h o \ 


(87) 


2n. 


(l-o 2 ) 


3/2 a. 


cos f 


(88) 


2n„ 


SN„ 


( X _e o 2 ) 


3/2 


Y 1_e o 2 sin f + 2 


GM 




r . T d r T 
— sin f -jI — cos f 


dg a 0 


(89) 
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and, as before, 


-o 2 a 


1 - e 


(! “ ®o 2 ) 


e n — sin f 


3/2 a 


with the new meaning of r and f. Setting 


M. 


(9 a 0 O" 


+ N. r 


(9 a Q £1" 
(97 


(i = 1, 2, 3) 


(90) 

(91) 


(92) 


and taking into consideration that the operators <9/<?g, 7 (<9/(97) are commutative with the operator 
K, we obtain the following differential equations for the contributions to Hansen's elements under 
the influence of the direct planetary effects: 


dAT 


+ oM 1 

<9 a 0 tr 

+ *>N, r 

cl a 0 $r 

dt 

rT i 


cl 7 

dAV 

rt 2 

+ bU 2 

* a o 

+ iN 2 r 

a 0 Q" 

dt 

<?g 

(97 


d h o 

dt A 17 " r T 3 - 

where the operator r is defined as 


(93) 

(94) 

(95) 


(96) 


and I is the identity operator. The values of n Q bz, and b\p taken from the solution of the main 
problem can be used in the computation of Equations 93 through 96. 

Our approach here differs from the approach to the main problem, in which we preferred the 
method of iteration to obtain the solution. In the main problem it is difficult to split the perturba- 
tions into different orders if the satellite theory is a purely numerical one. Generally speaking, 
the second iteration cycle only starts to produce the meaningful approximations to the amplitudes 
and the mean motions of the arguments of the periodic terms in the expansions of the coordinates. 

The uniformity of programming is an additional important factor which favors the use of the itera- 
tion process in the main problem. In the theory of the direct planetary effects the situation is dif- 
ferent, because (m"/l + m")(n 0 "/n 0 ) 2 9 the basic small parameter, is much smaller than the correspond- 
ing small parameter of the main problem. 
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The most significant part of the direct planetary effects is associated with the elliptic approx- 
imations to the motion of the sun and of Saturn. The terms which are linear with respect to per- 
turbations in the right sides of Equations 93 through 96 are the only terms which might produce 
some small but noticeable additional effects if the computation is performed with an accuracy of 
10~ 8 to 10~ 9 . In practice the operator r, like K, will always be reduced to only a few terms. The 
purely secular part in Equation 94 must be transferred to the corresponding equation in the main 
problem in order to determine the mean motion Jn 0 y dt of the pericenter relative to the ideal 
system of coordinates. 

THE DIRECT PLANETARY PERTURBATIONS OF THE POSITION 
OF THE ORBITAL PLANE OF THE SATELLITE 

Previously (Musen, 1959) we have established the following equations governing the oscilla- 
tions of the orbital plane around its mean position: 




1 h 

a o n o 

a(a 0 ft) 

dt ' 

+ n 0° X 2 

+ 2 ( x + ^ TT 0 

f 1 _e o 2 

dZ 

d^ 


1 h 

a o n o 

<9(a 0 ft) 

dt 

n o " -i 

4 2 ^ + ^) h 

- e o 2 

dZ 

d *3 


1 h 

_ a o n o 

f ' y ( a o s ‘) 

dt 

+ n 0 T ''4 

+ 2 + h 0 

t 1 " e o 2 

rlZ 



1 h 

a o n o 

<9(a 0 ft) 

dt 

— n t N 
1 1 0 ■ 3 

‘ 2 C 1 4 *0 h~ 

A IIq 

A^o 2 

dZ 


X 4 ^ COS ( f +r O ' ^3 i| Sin < f 


■-0 


“ v 3 cos (f + '*> - -4 i7 s i n ( f + '* ) 


] 


cos ( f * o ) 4 — sin(f-a) 

1 rf o 


’>[ ± - 7-1 

~ > - 1 a Q cos ( fil ) + S a 0 sin(f fa.) 


(97) 


(98) 


(99) 


(100) 


where aft :, 7. designates the component of the disturbing function normal to the orbital plane. In 
our present case, 


rift" /I 1 \ 

<9Z - m \| D -r 13 ~ D 3 / C 


( 101 ) 


or, in the expanded form, 


i? tt” 
clZ 


T m' C 


3r • D [iS (r • D) 2 

D 5 _ 2 D 7 


3 r 2 


D 5 J 



( 102 ) 


At this point it is convenient to change the meaning of notations in the same manner as we did in 
the previous section. After this change we can set 


r • D - A( r cos f) + B(rsinf) 
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and from Equation 102 we deduce with the new meaning of the notations, 


*( a >rr) 

SZ' 


+ 


3m" / n o \ / a o \ ^ AC 

l+m"\ n o/ \ ® / \a Q " 2 


r BC 7 . - 

— cos f + — — sin f 


We set 



(103) 


1 n 0 <9 ( a o n ") f 
A c - J —= jz — cos ( f + to) 

r 1 - e o 0 


(104) 



'' [ a oj - ' ) r 

r?Z 


— sin (f 4 'O 


(105) 


Then Equations 97 through 100 for the contributions to \ (i = 1, 2, 3, 4) as caused by the direct 
planetary effects, after some easy transformations, take the form 


cl * \ 


dt 



\ cos 


0 / 

2 ~ ( V 


■\. ~ \ 

4 s 


dAA 2 
dt" 



A cos 


- (A, 


s * + A 

3 s 


oK 


.) ■ 


dA \ 3 

dt 


r + v + 


IT - 1 


0 / \ 

A sin -y + (A ft \ A ft A , ) 

s z ' c 2 si/ 


dA\ 4 
~ dtT" 


- r + v + 



° / 

A c sin y - ( A ( 


^1 " A s 6A 2 ) • 


(106) 

(107) 

(108) 
(109) 


The purely secular terms should be transferred from Equations 107 and 108 to the corre- 
sponding equations of the main problem for the purpose of determination of the mean motion of 
the node and of the perigee. Normally the force component normal to the orbital plane is smaller 
than the components in the orbital plane. Thus the terms in the right sides of Equations 103 
through 109 which contain the perturbations as factors will be very small and in most cases can 
be omitted. 
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REQUIRED MODIFICATIONS IN THE MAIN PROBLEM 


In this section we discuss the changes in the main problem which are necessary to account for 
the deviation of the solar motion from the elliptic Keplerian motion. The expansions of the dis- 
turbing function and of its derivatives as given in the author T s previous work (Musen, 1963) require 
some modifications. The factor 1 + v' must now be attached to 7', and 7', f ' must now be con- 
sidered as functions of the disturbed mean anomaly V . Furthermore, to make the arguments 
linear with respect to time from the outset, we must expand the disturbing function and its deriva- 
tives in powers not only of n Q Sz as before, but also in powers of n ' Sz\ To perform this ex- 
pansion we must distinguish between the two uses of the symbols g and g' ; that is, between their 
use as arguments in the expansion of the perturbations and their use as the TT elliptic M g and g' 
(constituting the main parts of l and V respectively). The expansion in powers of n Q Sz, n Q ' Sz' 
requires formation of the derivatives with respect to the elliptic g and g' only. 

We resort to the standard Hansen device by introducing the temporary notations y, y‘ for g 
and g' respectively. After the expansion in powers of n Q Sz and n Q ' Sz' has been completed, we 
remove this distinction and return to the notations g and g' . We generalize here the authors 
previous expansion (Musen, 1963) by introducing the disturbing function 0* of the same external 
form as in the previous work but with g ,y and g' , separated. We have 


a o fi r 


a 0 °2 


a o ‘V 


where 



( 110 ) 


( 111 ) 


( 112 ) 


(113) 


The symbols s* , p*, q* are defined by 


* 

s 

+ (\ 2 -s 2 K* - * a * + iy 2 -y 2 ) s * - 2^ \ < ■ 

(114) 

* 

1 + V P a 0 


p 

1 + v‘ a o o' 

(115) 

* 

1 a o' 

(116) 

q 

1 + v’ p‘ 
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where 


s 


* 

1 


s 


* 

2 


s 


* 

3 


S 


* 

4 


1 + v p a o 
1 + V 1 a o p‘ 

1 + V p a o 
1 + l/ # a o p' 

1 + 1 / P 9 o 
1 + V ' a o P ' 

1 + 1/ p a 0 
1 + L/ ' a 0 p' 


COS (0 + (p‘ + CO T to* ) , 


sin (0 + 0 1 + Co + co' y , 


COS ( 0 ~ 0 ' + <*> - a/ ) , 


sin(0-0' + uj - oj ' ) , 


( 117 ) 


( 118 ) 

( 119 ) 


( 120 ) 


in which ,o, j, f ; , ct' are defined by the standard formulas 


/ ’ 

— cos 0 ~ cos t - e , 

0 0 


k sinrp 


j/Y ~ e 2 s i n € 


t - e Q sine = / , 


p' 

— - cos 0 - cos e ' - e ' , 

a o ° 


P' 



s i n 0 ' 



sine' 


e ' - e ' sin e ' - y ' . 


The expansions of 


p cos 

sin ^ ’ 


P 


a 


o 


and of 


a Q COS 
p' sin0, ’ 
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in terms of 7 and 7 ' respectively, can be obtained either by using Bessel functions or by means of 
the numerical harmonic analysis. After these expansions have been obtained, we can easily deduce 
the expansions of Equations 111 through 120 in terms of the arguments g, g' , g" , co 9 a/ . 

We now set 


where 


a 0 n = T' a 0 n* , 


e a 0 n 
p dp 


( n + 1) T' a Q SI* , 


T' 


L 


( n o' * z ')" _jT_ 

n! dy' n 


( 121 ) 


is the Taylor operator, and the "bar-operator" designates the replacement of 7 ' by g\ This 
modified value of the disturbing function must be used in association with the formulas and the 
differential equations developed for the main problem instead of its previous value as given in the 
earlier work (Musen, 1963). The form of the differential equations and of the integration procedure 
undergoes no change. However, the form of the solution will differ from the solution to the main 
problem. The series representing the coordinates will no longer be purely periodic. The inclusion 
of the perturbations of the sun and of Saturn will introduce very small mixed terms into the expan- 
sion of the coordinates. The mean motions of the node and of the perijove become the power series 
(in practice, polynomials) with respect to time. Consequently, Hansen's long-period arguments will 
become 

a- = o 0 4 n 0 (y + a p -n) dt , 

•'o 

" ' r 7 + n Q (a + 77 + y' ) dt . 

■'o 


CONCLUSION 

We have proposed here the theory and computational scheme which permit addition of the direct 
and the indirect planetary perturbations, with a high degree of accuracy, to the author's modifica- 
tion in 1963 of Hansen's theory. We suggest application of the method presented here to the deter- 
mination of the Saturnian effects in the motion of the outer Jovian satellites. No precise statement 
can be made in advance concerning the selection of terms. Each Jovian satellite displays its own 
peculiarities of motion, so the final selection of the significant periodic terms must be left to the 
electronic computer. The expansion of the disturbing function and of its derivatives, also the 
differential equations of the problem, are reduced to a form convenient for programming. 
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